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Abstract 

We have searched for the baryon-containing radiative penguin decays B~ — > Apj and B~ — ► 
S°p7, using a sample of 9.7 x 10 6 SB events collected at the T(4S I ) with the CLEO detector. We 
find no evidence for either, and set 90% confidence level upper limits of 

[B(B~ -> Ap 7 ) + 0.38(5" -> S°p 7 )] S7>2 . GcV < 3.3 x 10" 6 , 

_» s°p 7 ) + 0.4B(B~ -» Ap 7 )]i? 7 >2.o GcV < 6.4 x lO" 6 . 
From the latter, we estimate 

B(B — > X S 7,X S containing baryons) e 7 >2.o GcV < 3.8 x 10~ 5 . 

This limit implies upper limits on corrections to CLEO's recent measurement of branching fraction, 
mean photon energy, and variance in photon energy from b — > 57 that are less than half the 
combined statistical and systematic errors quoted on these quantities. 



*Present address: McGill University, Montreal, Quebec, Canada H3A 2T8 
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The branching fraction for the radiative penguin decay b — * 57 has been shown to place 
significant restrictions on physics beyond the Standard Model (SM) pj. The photon energy 
spectrum, in contrast, is insensitive to beyond-SM physics |2j, but provides information on 
the b quark mass and momentum within the B meson, information useful for determining 
the CKM matrix elements \ V n b\ and \V c b\. Some measurements of B{b — > sj) including 
the most precise one to date p, and the best measurement of the photon energy spectrum [5J] 
include a technique ( "pseudoreconstruction" ) that has reduced sensitivity to those B — > X s ^ 
decays with baryons in the final state. It is therefore important to determine what fraction 
of B — > Xg^ decays lead to baryons, or to place an upper limit on that fraction. 

Calculations of the mass distribution of the sq spect ator system that hadronizes into X s 
|^ show 1/3 of the spectrum above 2.05 GeV/c 2 , the threshold for Ap (the lightest 
baryon-containing final state), so a sizeable rate for b — > 57 with baryons would not be 
unexpected. Of the spectrum above Xp threshold, 2/3 is below 2.5 GeV/c 2 , so one expects 
the baryon-containing final states to be dominated by AN and EAT. Thus, measurements 
of the branching fractions for B — > ANj and B — > EA^y would help estimate a correction 
to the b —>■ sj branching fraction and photon energy spectrum. In addition, the decay 
B — > AN'-f provides a method for determining the helicity of the photon in b — > 57. For 
AA^ systems near threshold (s-wave), or for A and N near back-to-back to the photon (thus 
orbital angular momentum perpendicular to the photon), the A and 7 have the same helicity. 
A measurement of the A helicity, via its decay angle distribution, gives a measurement of 
the 7 helicity. We have therefore conducted searches for B~ — > Ajry and B~ — > E°p7 and 
their charge conjugates. 

The data used for this analysis were taken with the CLEO detector Q at the Cornell 
Electron Storage Ring (CESR), a symmetric e + e~ collider, and consist of 9.1 fb _1 on the 
T(45) resonance (9.7 x 10 6 BB events) and 4.4 fb" 1 60 MeV below the resonance. We select 
hadronic events that contain a A — ► pii~ , a p, and a high energy photon (E l ^ b >1.5 GeV), 
or contain a A, a p, and a high energy photon. (Henceforth, charge conjugate modes are 
implied.) For the B~ — > T,°p~f search we do not reconstruct the E°, but analyze as if the 
decay were B~ — > Apry, not detecting the soft photon from E° — > A7. The high energy 
photon must lie in the central region of the calorimeter (| cos# 7 | < 0.7), must not form a tt° 
or 77 meson with any other photon in the event, and must have a lateral energy distribution 
in the calorimeter consistent with that for a photon. The A requirements involve significance 
of displacement of vertex from interaction point and consistency of dE/dx and time of flight 
of the decay proton candidate with expectation. They result in a A candidate sample that 
is 90% pure. The antiproton candidate must pass dE/dx and time of flight requirements 
and must not form a A with any 7r + candidate in the event. 

We compute the standard B reconstruction variables M can( j = J -Eb eam — -P C and anc ^ = 
-Ecand — -^beam, keeping for further analysis events with M can( j > 5.0 GeV/c 2 and \AE\ < 
0.5 GeV. Pcand and E cand are computed from A, p, and 7 only, both for B~ — > Ajry and 
B~ — ► S°p7. With this event selection, there is negligible background from other B decay 
processes, but substantial background from continuum processes: initial state radiation, 
photons from decays of 7r° or rj that have escaped the veto, photons from decays of other 
hadrons. To suppress the continuum background, we compute twelve event shape variables 
and apply loose cuts on three of them. The twelve variables are then used as inputs to 
a neural net. The net is trained to distinguish between signal and continuum background 
using Monte Carlo samples of each. Monte Carlo samples distinct from those used to train 
the net are used to determine that cut on neural net output which would give the lowest 
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upper limit on the branching fraction should the branching fraction actually be zero, and 
also that cut which would allow us to see the smallest possible signal. Those two cuts differ 
little, and we use their average. 

The event shape variables are calculated in two frames of reference, the lab frame and 
the primed frame, the frame of the system recoiling against the photon. Variables in the 
primed frame are better at rejecting initial state radiation; those in the lab frame are better 
at rejecting other continuum events. The twelve input variables to the neural net are 1) 
| cos 9 a \, where 9 tt is the angle between the thrust axis of the candidate B and the thrust 
axis of the rest of the event, calculated in the lab frame; 2) | cos6^|, the same, but calculated 
in the primed frame; 3) the thrust of the candidate B; 4) the thrust of the rest of the event; 

5) i?2, the ratio of second and zeroth Fox- Wolfram (if moments, calculated in the lab frame; 

6) R' 2 , the same, but calculated in the primed frame; 7) | cos^'l, where 9' is the angle between 
the photon and the thrust axis of the rest of the event, all calculated in the primed frame; 
8) and 9) energies in 20° and 30° cones about the photon direction (excluding the photon 
energy); 10) the sum of the magnitudes of the component of momentum perpendicular to 
the thrust axis of the candidate B, for particles more than 45° from this axis, divided by 
the sum of the magnitudes of momentum for all particles, in both sums excluding particles 
from the candidate B, in the lab frame; 11) the same, but evaluated in the primed frame; 
and 12) cos9b, where 9b is the angle between the beam direction and the direction of the 
candidate B. 

The loose cuts are R 2 < 0.5, R' 2 < 0.3, | cos 6^1 < 0.8. Having obtained substantial 
suppression of background with the loose cuts and the cuts on the net output, our final 
selection is from the 2D distribution in M can( j — AE space. We define a "signal box" 
\AE\ < 84 MeV, |M cand - M B \ < 8 MeV/c 2 , which, based on Monte Carlo simulation, 
should contain ~90% of the B~ — > Apj signal events and (0.75 ± 0.15)% of the background 
events. We use the yield of events in the large M can( j — AE region (excluding the signal 
box), M can d > 5.0 GeV/c 2 , \AE\ < 0.5 GeV, to predict the background in the signal box. 
For B~ — > S°p7, we shift the signal box by 114 MeV to negative AE, compensating for 
the missing soft photon from S° — > A7. The shifted signal box should contain ~80% of the 
B~ —>■ S°p7 signal events. 

The 2D distributions in M can( j — AE space, On-4S resonance and Below-resonance, are 
shown in Fig. ^ There are 84 events On-resonance, and 43 events Below-resonance (with ~ 
half the luminosity), leading to a background prediction of 0.6 events On, 0.3 events Below, 
in either signal box. In the B~ — > Apy signal box, we observe zero events On and one event 
Below. In the B~ — > TPpy signal box, we observe one event On and zero events Below. The 
one On event has a B rest frame photon energy of 2.18 GeV, estimated by imposing the 
constraints that the undetected S° — > A7 decay photon brings AE to zero, and combines 
with the A to give the S° mass. Thus, we have no evidence for B~ — > Apy, and have a 90% 
confidence level upper limit on its true mean of 2.30 events. For B~ — > S°|?7, with one event 
observed and a background of 0.6 events expected, we also have no evidence of signal. We 
use the pre- Feldman- Cousins PDG procedure (§[ for calculating upper limits. We allow for a 
systematic error in background by conservatively using only half the expected background 
in the upper limit calculation. That gives a "conservative 90% confidence level" upper limit 
of 3.64 events. With the additional requirement that the B rest frame photon energy be 
greater than 2.0 GeV, the background in the large M can d — AE region drops to 27 events 
On and 15 events Below, with 0.21 background events predicted for the B~ — > TPpy signal 
box. This leads to an upper limit of 3.80 events for E^ > 2.0 GeV. 



4 



*1 850403-001 



; — t* — | — i — | — r- 

Below-Resonance Data 



0.4 



0.2 



> 

2. o 

LU 
< 



-0.2 



-0.4 



On-Resonance Data* 



± 



5.00 5.05 5.10 5.15 5.20 5.25 5.00 5.05 5.10 5.15 5.20 5.25 

M B (GeV/c 2 ) 



FIG. 1: Mb — AE for On and Below resonance data and photons with El?° > 1.5 GeV. The 
solid box shows the tight signal box (5.272 < M B < 5.288 GeV/c 2 , \AE\ < 0.084 GeV) used for 
determining the B~ — » A7J7 yield. The dashed box is shifted downward in AE by 114 MeV and is 
used for determining the B~ — > S°p7 yield. 



The upper limit on the branching fraction will be those upper limits on the number of 
signal events, divided by the number of charged B's, by the A — > pir~ branching fraction, 
and by the detection efficiency. 

We assume equal number of charged and neutral B's, noting that a correction for this 
assumption can be applied at such time as the B + B~ to B°B° ratio in T(4S f ) decays has 
been well determined. Thus, we have 9.7 million charged B's. We assign a ±2% uncertainty 
to that number. 

We use Monte Carlo simulation to determine the efficiency for detecting B~ — > Ap~f, and 
B~ — > £°£>7. The Standard Model predicts left-handed photons and s quarks in b — > 57, 
and thus the A's in B~ — > Apy will tend to be left-handed. This tendency decreases the 
number of pions from A — > pn~ that decay against the A boost direction. Such pions are 
soft and more difficult to detect. Thus, SM decays will have a detection efficiency higher 
than that of unpolarized A's. For an upper limit, we conservatively assume unpolarized A's. 
For B — > S°P7, we also assume unpolarized A's. 

The efficiency as a function of Ap mass is shown in Fig. [21 The sharp fall-off near 3.5 
GeV/c 2 is caused by the photon energy requirement, E l ^ b > 1.5 GeV. The gentler decrease 
from 2.4 to 3.4 GeV/c 2 is caused by the background suppression requirements. Similar 
results are obtained for We assume a Ap mass distribution mass distribution) 

given by the parton- level hadronic mass distribution 0, E| times a phase space factor P/M, 
P being the momentum of A or p (E° or p) in the Ap rest frame, for that value of Ap 

(£°p) mass M. We have also used a weighting P 3 /M, appropriate for a p-wave system. 
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FIG. 2: Efficiency as a function of M A p, for signal Monte Carlo, with all cuts applied. 



As we are primarily interested in decays with a high energy photon, we compute the 
efficiency for the subset of events with B rest frame photon energy E 1 > 1.5 GeV (M\ p < 3.5 
GeV/c 2 ), and with £ 7 > 2.0 GeV (M Ap < 2.6 GeV/c 2 ). For B~ -> Ap 7 , for events with 
E 1 > 1.5 GeV we find efficiencies of 11.6% (for P/M) and 10.5% (for P 3 /M); for events 
with Ej > 2.0 GeV we find an efficiency of 12.4% in both cases. For the E 1 > 1.5 GeV case, 
we conservatively use the smaller efficiency. For B — > S°P7, for events with E y > 1.5 GeV 
we find efficiencies of 9.4% (for P/M) and 8.2% (for P 3 /M); for events with E 1 > 2.0 GeV 
we find an efficiency of 10.6% in both cases. 

There are also systematic errors in the efficiency from uncertainty in the simulation of 
the detector performance (track-finding, photon-finding, vertex-finding, resolutions) and an 
uncertainty in the modeling of the other B. We estimate these at ±8.2%. 

We obtain a conservative 90% confidence level upper limit on the branching fraction by 
using unpolarized A's, using the P 3 /M option for the Ap (S°p) mass distribution, and then 
increasing the limit so obtained by 1.28 times the quadratic sum of the two remaining sys- 
tematic errors, ±2% from number of B's and ±8.2% from detector simulation and modeling 
of the other B. 

While our specific goal in the first search was B~ — > Ap7, we also have sensitivity to the 
decay B~ — > S°p7, S° — > A7 in that analysis. Our efficiency for the latter decay is 0.3 times 
that of the former. Similarly, while our specific goal in the second search was B — > S°p7, 
we also have sensitivity to B~ — > Ajry, 0.4 that for B~ — > TPp^y. Hence, our primary results 
can be written 

[B(B~ -> Ap 7 ) + 0.3B(B~ -> £°p 7 )b 7 >i.5 GeV < 3.9 x lO" 6 , 



[B(B~ -> Ap 7 ) ± 0.3B(B- -> S°p 7 )] S7>2 . GcV < 3.3 x 10^ , 
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[B{B- -> £°p 7 ) + 0.4B(S" -> AP7)]s 7 >i.5 Gev < 7.9 x lO" 6 , 

-> £°p 7 ) + 0.4B(S" -> Ap7)]s 7 >2.o Gev < 6.4 x 1(T 6 . 

From these upper limits, we would like to obtain an upper limit on the branching fraction 
for b — > S7 leading to baryons. Our first step in this direction uses isospin considerations. 
The parton- level final states, su and sd, form an isospin doublet, and the hadronization 
process should conserve isospin. This gives B[B~ — > Apj) = B(B° — > Anj), and B[B~ — > 
£°p 7 ) = B(B° -> £°n 7 ) = 1/2B(B~ -> E'n-f) = 1/2B(B° -> S + p 7 ). Thus -> 
EN 7 ) = 3B(S" -> S°p7), and S(£ -> (A or E)iV 7 ) = 3B(B" -> S°p 7 ) + -6(5- -> Apy). 
Multiplying the last upper limit given above by 3, we have [B(B —>■ SA^ 7 ) + \.1B(B — > 
AA^ 7 )]£; 7>2 .o Gev < 1-9 x 10 -5 , which we use as our limit on B(B — > (A or S)iV 7 ). 

The above branching fraction limit is for both baryon and antibaryon in the lowest lying 
baryon SU(3) octet. We must also consider decays with one of the baryons in the decuplet 
(i.e., B — > SA 7 and B — > E(1385)iV 7 ), decays involving higher-mass octet and decuplet 
members, and non-resonant decays such as B —>■ (A or S)(iV or A)7r 7 . The requirement 
that Ej(B rest frame) be greater than 2.0 GeV translates into an upper limit on the mass 
of the baryon-antibaryon system of 2.60 GeV/c 2 . The various mass thresholds are: AN, 
2.05 GeV/c 2 ; EJV, 2.13 GeV/c 2 ; EA, 2.43 GeV/c 2 ; E(1385)iV, 2.32 GeV/c 2 . Thus, phase 
space will suppress the octet-decuplet rates relative to the octet-octet rates. Combining 
this with the falling parton-level hadronic mass distribution given by the spectator model 
0, or the calculation of Kagan and Neubert 0, we estimate a suppression of a factor of 
~4. This is partially compensated by the factor of 2 more spin states available in the octet- 
decuplet combination. A plausible assumption is that the octet-decuplet contribution would 
be ~l/2 that of the octet-octet contribution. Octet-decuplet pairs with an excited member 
are above the 2.60 GeV/c 2 cutoff imposed by the 2.0 GeV photon energy requirement, 
and octet-octet pairs with excited members have thresholds very close to the cutoff. Non- 
resonant (A or S)(iV or A)tt states will have thresholds below the cutoff, but will be phase- 
space-suppressed relative to (A or E)iV. From all this, we take as our working assumption 
B{b — > sj, with baryons) El >2.o GeV = 2B{B — > (A or £)iV 7 )£ 7> 2.o Gev, and hence B(b — > 
sj, with baryons) El >2.o Gey < 3.8 x 10" 5 . 

CLEO's recent study 5] of b —> 57 reported a branching fraction for E 1 > 2.0 GeV, 
corrected for the b — > dj contribution, of (2.94 ± 0.39 ± 0.25) x 10~ 4 . Our upper limit on 
the branching fraction for b —>■ sj leading to baryons, with E 1 > 2.0 GeV, 3.8 x 10~ 5 , is 
13% of that number. The recent study [5J] had an efficiency for detecting B — ► baryons 7 
that was 1/2 of that for modes not involving baryons. That implies an upper limit on 
the correction needed for the branching fraction reported there of 6.5%, less than half the 
combined reported statistical (±13%) and systematic (±8%) errors. 

CLEO's recent study [^] of b — > 37 also reported information on the photon energy 
spectrum: an average energy (E^) = (2.346 ± 0.032 ± 0.011) GeV, and a variance ((E y — 
(£ 7 )) 2 ) = (0.0226±0.0066±0.0020) GeV 2 . Both averages were taken only for photons above 
2.0 GeV. The average energy of B rest frame photons from events with baryons (averaging 
only for photons above 2.0 GeV) is ~2.1 GeV, 250 MeV lower than the published mean. 
The upper limit on the correction to the first moment is thus 6.5% of 250 MeV, i.e. 16 MeV 
(compared with the published statistical and systematic errors of 32 MeV and 11 MeV, 
respectively). The limit on the correction to the variance is 0.0025 GeV 2 , which is 36% of 
the combined quoted statistical and systematic errors on the variance. 
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In conclusion, we have conducted searches for the exclusive radiative penguin decays 
B~ — > Ap7, and B~ — > S°p7, found no evidence for either, and placed upper limits on them 
of 
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-6 



With plausible assumptions, this leads to the conclusion that b — > 57 decays with baryons in 
the final state and E y > 2.0 GeV constitute at most 13% of all b — > decays with E n > 2.0 
GeV. With that limit, the upper limit on corrections to our recent measurement [^] of the 
b — > 57 branching fraction, the mean energy of the photon, and the variance in the photon 
energy, are less than half of the combined quoted statistical and systematic errors. 
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